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FREQUENCY DOMAIN ECHO CANCELLER 

FIELD OF THE INVENTION 

The invention relates to digital communication systems. One application of 
the invention is systems that use multicarrier modulation methods, such as DMT (Discrete 
Multitone ) where DSL (Digital Subscriber Lines) twisted pair subscriber Imes are 
employed for high speed data transmission. Another application of the invention is 
OFDM (Orthogonal Frequency Division Multiplexiag) based systems. 

BACKGROUND AND SUMMARY OF THE INVENTION 

Referring to Fig. I, a DSL modem consists of a transmitter 12 and a 
receiver 14 connected to the subscriber line's twisted-pair 18 (telephone line) via a so- 
called hybrid circuit 16. The hybrid 16 is a two-toTour wire interface that connects the 
twisted-pair subscriber liae 18 to the four-wire modem (two-wire transmitter and two-wire 
receiver). An ideal hybrid is exactly matched to the impedance of the subscriber line and 
transports the transmitter signal only to the subscriber line and the received signal only to 
the receiver. 

A problem with matched hybrids (non-adaptive) is that different subscriber 
hues have different impedances^ and temperature variations make the impedances slowly 
time-varying. Consequently, the hybrid cannot be exactly matched to the impedance of 
the subscriber line, and the transmitted signal "leaks" back into the receiver through the 
hybrid creating an echo signal as shown in Fig. 1. From the receiver's point of view, the 
echo signal acts like a powerful noise source reducing the signal-to-noise ratio (SNR) and 
therefore reduces the performance of the receiver. An echo canceller may be used to 
reduce the influence of the echo in the received signal. 

The main purpose of the echo canceller is to reduce the echo signal 
produced when the transmitted signal "leaks" back into the receiver via the hybrid. The 
echo canceller works by first estimating the echo signal from the transmitted signal, and 



then subtracting the estimated echo signal from the received signal. In order to estimate 
the echo signal, the echo canceller must model the echo path channel. Generally speaking, 
the echo canceller performs a ''system estimation/' where the system to be estimated is the 
echo path channel. Due to the slowly time-varying nature of the echo path channel, the 
5 echo canceller must continuously, but infrequently, update the echo path channel 

estimation. This update is typically accomplished with an adaptive finite impulse response 
(FIR) filter to continuously track the slowly time- varying echo path channel. 

The echo canceller operates in two different modes: the training mode and 
the echo cancellation mode. In the training mode, the far-end transmitter is forced to be 

10 silent while the near-end modem transmits a training sequence in order to estimate the 

echo path channel. The estimation of the echo path channel may be accompHshed using a 
fast training algorithm followed by an adaptive algorithm such as the LMS-algorithm. 
After the echo canceller finishes the training mode, it jumps iato the echo cancellation 
mode where the echo canceller performs echo cancellation, i.e. it subtracts the echo signal 

15 from the received signal. 

Discrete Multi-tone (DMT) modulation is one form of multicarrier 
modulation that is attractive because of its ability to utili2e, and be implemented by, 
efficient digital signal processing techniques. The basic idea of DMT modulation is to 
partition the communication channel iato orthogonal and memoryless subchannels using 
20 the digital signal processing technique known as the Fast Fourier Transform and the 

Inverse Fast Fourier Transform (FFT/IFFT). Data transmitted through one subchannel is 
independent of other subchannels under some orthogonal conditions that will be 
described later. DMT also allows for efficient and flexible allocation of power and bits to 
each subchannel. 

25 The DMT modulated signal can be described as a composition of N/ 2 

independent carriers, where each carrier has equal bandwidth and a center frequency /, 
where / = / . . .N/2. For example, ADSL uses 255 carriers (subchannels) downstream and 
31 carriers upstream. Each one of the carriers is a Quadrature Amplitude Modulated 
(QAM) signal with an assigned amplitude and phase. The assigned values for the 
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amplitude and the phase are determined by a signal constellation containing a set of 
complex values. During start-up, the DMT-system measures the signal-to-noise ratio in 
each subchannel and assigns different numbers of bits to each subchannel (earner) to 
maximize performance. This process is known as adaptive bidoading. A subchannel with 
5 low signal-to-noise (SNR) is assigned a small number of bits or no bits, and a subchannel 
with high SNR is assigned many bits. This technique is robust in a typical Digital 
Subscriber Line scenario where the line conditions are unknown and slowly time-varying. 

A DMT signal diat is filtered through a physical channel will be subjected to 
inter-symbol-interference (ISI) because all practical physical channels have a memory. 

10 (The impulse response of the channel is non-zero for more than one point in time). This 
impairment of the physical channel causes the transmitted DMT symbols to interfere with 
each other (ISI). Furthermore, the transitions between the DMT symbol cause transients 
in the received signal. These transients cause mterference between the carriers in the same 
DMT symbol. This type of interference is known as inter-channel interference (ICI), 

15 where the "channel" refers to the subchannel. The ISI and ICI imply that the received 
subchannels do not remain orthogonal (mdependent) to each other, and therefore, the 
received data wiU be corrupted. 

The received echo signal, like the received data signal, will also be affected 
by the ISI and ICI, because all practical physical echo padi channels have memory. Thus, 
20 the echo from one carrier will leak into every other carrier within the same DMT symbol 
and the next transmitted DMT symbol. 

One way to circumvent, or at least reduce, the ISI and ICI in a DMT-based 
system is to add a cyclic prefix (guard time) to the beginning of each transmitted DMT 
symbol. The cyclic prefix (CP) is a mechanism to make the subchannels independent of 
25 each other, or m other words, memory-less. A cyclic prefix of L samples means that the 
last L samples of the N samples long time domain symbol is copied to the beginning of 
the time domam symbol. Thus, the total length of the prefixed time domain symbol is 
L + N samples. Figure 3 illustrates the cyclic prefixing. 



If the length of the cyclic prefix is greater than the length of the channel 
memory, then aQ ISI and all transients causing ICI will be contained within the cyclic 
prefix. The receiver can then remove all ISI and ICI by simply discarding the cyclic prefix. 
This means that the channel will only affect the ampHtude and the phase of each 
transmitted carrier, and different subchannels will remain independent (no ISI or ICI). 
Unfortunately, increasing the length of the cyclic prefix reduces the data rate. For 
example, suppose that one time domain symbol is N samples long, and the length of the 
cycUc prefix is L samples. The data rate is reduced by a factor of N/{N + L). 

In a practical DMT-system the echo path channel is much longer than the 
cyclic prefix, and therefore, the received echo signal is subjected to both ISI and ICI. 
Hence, the part of the echo signal that corresponds to both ISI and ICI must be estimated 
and subtracted as well. The present mvention provides an echo canceller that takes into 
account intersymbol interference and interchannel mterference when estimating an echo 
signal to be removed from a received signal. The echo canceller calculates the received 
echo signal completely in the frequency domam, which is particularly beneficial in a DMT 
type system because the transmitted data is akeady available in the frequency domain. 
Accordragly, the echo estimate is determiaed using a frequency domaia model of an echo 
path channel that includes the effects of interference like ISI and ICI. 

In a first example embodiment, the frequency domain model of the echo 
path channel is determined using a first matrix of coefficients and a second matrix of 
coefficients. The first matrix is combined with a current transmitted symbol, and the 
second matrix is combined with a previously transmitted symbol. The sum of these two 
combinations is used to estimate the echo signal. The coefficients of the first matnx 
represent how an echo from a currently transmitted frequency domain signal affects a 
received signal. The coefficients of the second matrix represent how an echo from a 
previously transmitted domain symbol affects the received symbol The coefficients of the 
first and second matrices are adjusted using a difference between the received signal and 
the estimated signal- 
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In a second example embodiment, the current transmitted symbol and the 
previously transmitted symbol are partitioned into real and miaginary parts before bemg 
combined with matrices as described in the first example embodiment- This operation 
reduces the computational complexity by a factor of two compared to the first example 
5 embodiment. 

In a third example embodiment, the currently transmitted symbol is 
combined with a first column vector. The previously transmitted symbol is multiplied with 
a complex exponential term, in order to compensate for the cyclic prefix, and then 
subtracted from the current transmitted symbol- The resulting signal is combined with a 
10 matrix. The vector and matrix combinations are summed and used to estimate the echo 
signal to be removed from the received symbol This third example embodiment reduces 
the computational complexity approximately by a factor of two compared to the first 
example embodiment. 

In a fourth example embodiment, the operations described in the second 
15 and third example embodiments are combined resulting in a reduction of computational 
complexity approximately by a factor of four compared to the first example embodiment. 

In a fifth example embodiment, when the transmitter of a transceiver has a 
lower sampling rate than its receiver, the received echo signal is interpolated at the 
receiver. The interpolation may also be combined with the operations described m the 
20 second example embodiment to reduce the computational complexity by a factor of four 
compared to the first example embodiment. 

In a sixth example embodiment, when the transmitter has a higher rate than 
the receiver, the echo signal is decimated at the receiver. The decimation may also be 
combined with the operations described m the second example embodiment to reduce the 
25 computational complexity by a factor of four compared to the first example embodiment. 

In a seventh example embodiment, when the transmitted symbols are not 
aligned in time with the received symbols or frames, an asynchronous echo canceller may 
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be used. The asynchronous echo canceller may be combined with any of the previously- 
described example embodiments. 

In a DMT type modem, because all of the data is already available in the 
frequency domain, there is no need to do any additional Fourier transform or mverse 
5 Fourier transform operations, except for the asynchronous echo canceller embodiment 
which requires one extra IDFT. Moreover, some of the coefficients in the matrices may 
be very small and ignored, further reducing the computational complexity and memory 
required. 

RRTKF DESCRI PTION OF THK DRAWINGS 

10 The foregoing and other objects, features, and advantages of the invention 

will be apparent from the foUowing description of preferred, non4imiting example 
embodiments, as weU as illustrated in the accompanying drawmgs. The drawings are not 
necessarily to scale, emphasis instead being placed upon illustratmg the principles of the 
invention. 

15 Fig. 1 illustrates a symmetric-rate DMT transceiver where the hybrid is 

unmatched to die subscriber line; 

Fig. 2 illustrates cyclic prefixing; 

Fig. 3 shows a symmetric-rate DMT transceiver with an echo canceller; 

Fig. 4 shows a symmetric-rate DMT transceiver with one example 
20 embodiment of a frequency domain echo canceller; 

Fig. 5 shows a flowchart of procedures performed by the echo canceller in 

Fig. 4; 

Fig. 6 shows a symmetric-rate DMT transceiver with another example 
embodiment of a domain echo canceller; 
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Fig. 7 shows a flowchart of procedures performed by the echo canceller of 

Fig. 6; 

Fig. 8 shows another example embodiment of the invention as applied to 
asymmetric data rate transceivitig situations; 

5 Fig. 9 shows a symmetric-rate DMT transceiver with another example 

embodiment of the invention as applied to asynchronous data transceiving situations; and 

Figs. 10 and 11 show the learning curve of the echo canceller of Fig. 6 
implemented m a DMT-based, symmetric-rate DSL transceiver, using a simple lowpass 
model of the echo path channel. 

10 DETAILED DESCRIPTION OF THE DRAWINGS 

In the following description^ for purposes of explanation and not limitation^ 
specific details are set forth, such as particular embodiments, protocols, data stmctures, 
and techniques, in order to provide a thorough understanding of the present iavention. 
However, it will be apparent to one skilled in the art that the present invention may be 

15 practiced in other embodiments that depart from these specific details. In other instances^ 
detailed descriptions of well-known methods, systems, and devices are omitted so as not to 
obscure the description of the present iavention with unnecessary detail. Moreover, 
individual function blocks are shown in some of the figures. Those skilled in the art wiU 
appreciate that the functions may be implemented using individual hardware circuits, using 

20 software functioning in conjunction with a suitably programmed digital microprocessor or 
general purpose computer, using an Application Specific Integrated Circuit (ASIC), and/or 
using one or more Digital signal Processors (DSPs). 

A discrete multi-tone modem 20 in which the present invention may be 
employed is illustrated in Fig. 3. The DMT modem includes a DMT transmitter 22 and a 
25 DMT receiver 26 coupled to the subscriber Hne via a hybrid 24. Serial input data are 
grouped into blocks, converted into parallel form^ and encoded by an encoder 28. Each 
one of N/ 2 active subchannels contains a number of bits allocated to that subchannel. 



8 



Each block of bits is mapped into a two-dimensional signal constellation point (i.e., a 
complex number). The output from the encoder 28 is N/2 complex numbers (ie,^ 
subsymbols), one for each block of bits. The N/2 complex subsymbols constitute a DMT 
symbol in the frequency domain and are extended to Hermitian symmetry (i.e.^ conjugate 
5 symmetry) before being input to block 30^ which calculates the inverse discrete Fourier 
transform (TDFT). Due to the Hermitian symmetry property^ the N-IDFT output is a 
sequence consisting of iV real samples. The output sequence consists of the superposing 
of the N/2 modulated orthogonal carriers. The parallel N-IDFT output is converted to a 
serial stream by a parallel-to-serial converter 32. The digital modulated data stream is 
10 cychcally prefixed in the Add CP block 34^ converted to analog form by digital-to-analog 
converter 36, filtered at the transmit filter 38, and passed to the hybrid 24 for transmission 
over the subscriber line. 

At the receiver end, the received analog signal from the hybrid 24 is filtered 
in receive filter 40, converted to digital form in the analog to digital converter 42, and 

15 stripped of the cychc prefix in the Remove CP block 44. The resulting digital signals are 
converted to parallel format in the serial- to-parallel converter 46, The parallel digital 
signals are demodulated in the N-DFT block 48, The frequency domain echo canceller 50 
estimates the echo caused by the reflected signals from the DMT transmitter 22, and those 
echo signal estimates are removed from the demodulated signal at block 52. The 

20 remaining frequency domain signal is equalized, decoded, and converted back to a serial 
form by decoder 54. 

The remaining description will focus on the frequenq^ domain echo 
canceller 50. To imderstand the mathematical foundation for the echo canceller of the 
present invention, a mathematical expression of the received echo signal in the frequency 
25 domain is derived as a function of the transmitted frequency domain signal and the echo 
path channel. A significant feature of the invention echo canceller is that the received 
echo signal is cancelled completely in the frequency domain. This is particularly beneficial 
in a DMT transceiver where the transmitted data is already available in the frequency 
domain. There is no need to do any extra DFT or ID FT calculations except for the later- 
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described asynchronous echo cancellation embodiment which requires only one extra 
IDFT. 

As explained above^ the received echo signal includes ISI and ICI because of 
the long impulse response of the echo path channel. In order to estimate the received 
5 echo signal^ a frequency domaia model of the echo path channel should include the effects 
of the ISI and the ICL 

The following notations are introduced: 

x, = [x,(o)iX,(i),...,x,(A^-i)r 0) 

denotes the i:th transmitted frequency domain symbol containing the JV transmitted 
10 frequency subsymbols ^where subscript i denotes the DMT symbol number and k 

denotes the frequency sample mmiber. Note that the vector Xi is Hermitian symmetric in 
order to obtain a real-valued sequence after the IDFT. 

The Hermitian symmetry of the vector Xi is defined as 

XX>t)-X;(iV-^),for^-A^/2 + l,A^/2 + 2,...,A^-l, (2) 

15 where * denotes complex conjugate. Furthermore, the DC term XfO) and the Nyquist 

term X (Njl) must both be real-valued. Note that X t is the transceiver's own transmitted 
signal that generates the echo signal Y?. The vector 

Y,=[>^(0U(l),...,F,(iV-l)r (3) 

denotes the i:th received frequency domain (Hermitian symmetric) echo symbol contammg 
20 the received frequency subsymbols Yi(k). OrJy the received echo signal is considered, and 
the far-end data signal is ignored. Thus^ the scenario emulates the echo canceller training 
mode where the far-end modem is silent. 

The iV~point IDFT of X i can be expressed by 

s, =QX„ (4) 
25 where Q is the IDFT matrix defined by 

r 1 1 i — 

[q]„.^=— ,for« = 0,l,...,7V-landA: = 0,l,...,7^-l. (5) 
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Defining as a matrix containing the last jL rows of Q, i.e. for n — 0^ 1^ 
. .jL -/ resxilts in the following: 

[QcX,=[Q]N-L.n,, (6) 

The transmitted, prefixed time domain symbol is written as: 



Q 



(7) 



The sequence consisting of the / -7:th and /.th prefixed time domain DMT- 
symbol can be written as: 







Q 










Q 



X, 



(8) 



where (Z)^ denotes a matrix with L rows and N columns consisting of zero-valued 
10 elements. 

The length of the echo path channel is assumed to be M samples long, and 
the impulse response of the echo path channel is denoted with die convolution matrix, G, 
which is of size (2N^2L) x (2N~h2L). 
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(9) 
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15 The received time domain echo signal equals 



(10) 
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In the receiver, only the i:th DMT-symbol is of iaterest. It contains both ISI and ICI. The 
cyclic prefix (CP) is discarded, and the received vector y is multiplied with the DFT matrix 
R of size (NxN), defined as 



-7 — nk 

N 



(11) 



With Rtor defined as R,„, = \0nx{n+l) (Dnxl rJ the received /:th frequency domain 
DMT-symbol (vector of size NxV) is obtained as: 



Y, - R,,GQ^,X - AX - [Ap,,, ] 



X 



(12) 



This can also be written as 

~ ^ prev J-1 current 

where the following matrices of size (N xN) are defined: 



(13) 



Po.o Po,i 
Pi.o Pu 



Po. 



Pn~,. 



(14) 



"0,0 ^0,1 



(15) 



Thus, the k:xh received frequency echo sample ia the /.th symbol can be 



written as 



N~l 



(16) 



W=0 
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From equation (16)^ the first term c^^X^(A:) describes the direct echo firom 

AT-l 

the carrier number k in the symbol number i. The second term ^^^k^m^X^ describes 



m=0 



the ICI horn, symbol number i The third term ^/^^^^.^^ (^) describes the combined ICI 

and ISI from the previously transmitted symbol mmiber i-l. If the cyclic prefix was longer 
5 than the echo channel memory, then only the first term would be non-zero, and there 
would not be ISI or ICI terms because the matrix A«,^^«/ will only have non-zero values at 
the main diagonal (upper left comer to the lower right comer), and the matrix A^;r^ would 
contain only zeros. 

Hence, the matrix Kprev describes how the echo from the previously 
10 transmitted frequency domain symbol leaks into or distorts the received symbol The 
matrix describes how the echo from the currently transmitted frequency domain 

symbol leaks into or distorts the received symbol- Accordingly, the two matrices A^w and 
Kmni model how ISI and ICI will affect the echo signal completely in the frequency 
domain. By determining the coefficients for the matrices Kprev and Kcumnt^ the echo in the 
15 frequency domain wiU be known and can then be subtracted from the received signal. The 
coefficients of the unknown matrices Aprev and Aairrent can be adaptively identified in order 
to estimate the received echo signal. Matrices Kprev and Kcmrent can also be calculated firom 
equation (12) by estimating the echo path channel impulse response (matrnc G). 

Fig. 4 shows a transceiver using Discrete Multitone (DMT) as a modulation 
20 scheme with a first, non-limiting, example embodiment of a matrix-based frequency 

domain echo canceller (MBAEC) stmcture. There are two complex matrices H/ and W, , 
numbered 62 and 64, respectively. Each matrix has a size of (N xN)^ where JVis the size 
of the IDFT and DFT. The received echo signal is estimated in the frequency domain by 
multiplying the matrix with the currendy transmitted frequency domain symbol (an 
2,5 N X 1 vector), X^, output by the encoder, and multiplying the matrix Wi with the previously 
transmitted frequency domain symbol (an N x 1 vector), X^.y, ou tput from delay 60. The 
two products H.X^and ^ are added together, i.e.. 
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(17) 



where 



H, = 



^0,0 ^0,1 
^1.0 ^1.1 



H 



N-\,N-\ 



(18) 



and 



0,0 0,1 
^,0 W,, 



w 



w w 



(19) 



Thus^ the k,th estimated frequency echo sample in the symbol is 



N'-l 



N-1 



t {k) = H,^,X, {k)+j:H,,„X, {m) + (m) (20) 

The hat symbol of ff,{k) indicates that it is an estimation of the received echo signal (k). 

The estimated echo signal is subtracted from the received frequency domain 
10 symbol, and the resulting error vector = Y ^ - □ ,is used by an adaptive algorithm^ e.g.. 
Least Mean Square (LMS), to adjust the matrix elements of H/ and Wi. In regular 
operation, (after training) the error signal is also the desired data signal. Adjusting and 
with the LMS-algorithm yields: 

H,,, =H,+//E,(x:/ (21) 

15 W,,,-W,+//E,(x:J (22) 

where - - = [E^ {o\ (l), ,.,EXN ~ l)Y is the i:th received frequency domain error 
vector and// is the step size of the LMS algorithm. One can prove that the LMS 
algorithm will adapt the coefficients of Mi to an estimate of A^^^^ and the coefficients of 
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W? to an estimate of A^^^^ . Most of the coefficients of the matrices and are 

complex-valued. Hence, all arithmetic computations are complex additions, subtractions, 
and multipUcations. 

Some of the coefficients (elements) in A^^^^^^ and Ap^.^^ will be almost zero 

and can therefore be neglected to reduce the computational complexity in equations (21) 
and (22). This wiU also reduce the amount of memory. In fact, Acurremand Aprevmay often 
be band matrices meaning that maialy the strong diagonal (left to right) and some of the 
upper and lower diagonals have non-zero values. The nimiber of diagonals to be taken 
into account depends on the echo path channel (ISI and ICI) atid the desired level of echo 
cancellation. 

The operation of an echo canceller in accordance with the first, non-limiting, 
example embodiment and illustrated m Fig. 4 is now described in conjunction with the 
MBAEC routine (block 70) shown in flowchart format in Fig. 5. A signal to be 
transmitted Xi is detected in the frequency domain (block 72). That signal is multiplied by 
matrix 62 and the previously transmitted frequency domain signal Xi-i is miiltiphed by 
matrix Wi (block 74). The matrix products are summed together, i.e., HiXi + WiXi-i, to 
estimate the current echo (block 76). The estimated echo is subtracted from the 

actually received symbol Yi to determine an error Ei (block 78). The coefficient values of 
the H and W matrices are then adjusted for the next signal period / + 1 with the detected 
error Ei using for example the LMS algorithm (block 78). 

When the transmitted time domain signal is real valued, e.g., in the case of 
DMT, it is possible to avoid the use of complex arithmetic. In order to obtain a real- 
valued, time domain signal, the frequency data must be Hermitian symmetric, and the DC 
and Nyquist terms must both be real-valued. In a second, non-limiting, example 
embodiment of the present invention, a real-valued implementation of the MBAEC 
reduces computational complexity by a factor of two compared to the general case 
described above. A real-valued vector is introduced as 
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(23) 



where the vectors X,,^,, and X,„^^ , are defined as 



^real , = [M^X Xreal. ,(1), Xreal, ,(2),. . ., Xreal, ,{N I 2 - 1), X,{N I 2)f 



(24) 



and 



^ima^ 1 ~ til), Ximag, i{2\,,., Ximag, i(N / 2 ~ 1)] 



(25) 



In equations (24) and (25), Xreui,,in) ^ Re{X(/2)}and X,,.ag,,(n) = Im{X:(/2)} are defined. 
Furthermore, X(0) is the DC term, and X(N/2) is the Nyquist term of X, . The vector 
X, contains exactly the same information as the vector X, dejSned in equation (1), but in a 
slighdy different order. 

The matrices Q,„, and Q,,^ are defined for ^ = 0, 1, . . . , N/2 - 1 and 
«^0,1,...,A''-1 as 



[Q^oJ t = — cos 

cos Jn,k 



kn 



, for k = 0X-;N/2 



(26) 



[Qs.L-|sin[^(A: + l)«j, fork^0X...,N/2-2 (27) 



Furthermore^ the matrix Q is defined as 



Q cos ^ sm 



(28) 



Equation (4) may be rewritten 



(29) 



Qcp is farther defined as a matrix containing the last L rows of Q , i.e. for « = 0, 1, . . . Z, - 1 
we have 
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[q4.^=[qL.«,. (30) 

The transmitted prefixed time domaia symbol may now be written as 



Q 



(31) 



As in equation (8)^ the sequence of the i-/:th and i:th prefixed time domain DMT symbols 
can be written as 



X = 







Q 






Q 


®NxN 


Q 



X, 



(32) 



With G defined in equation (9), the received time domain signal may be written as 

y-GQ,,,X (33) 

In the receiver^ we are only interested in the /:th DMT symbol which 
contains both ISI and ICI. The cyclic prefix is discarded, the received vector y is 
multiplied with the matrix R^^^ defined as 



^^^^ ^[0Nx(n^l) QnxL r] 



(34) 



where 



R-[Qcos -Qsm] (35) 

The received /:th firequency domain data vector is now expressed as 

Y, = R,„,GQ,,X = AX = ] . (36) 



X. 



where we have defined 
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Y = 



Y 

real, J 



, 7mag, I 



(37) 



and 

^reaU = '0), ^real, ,(2),..., Fw, ,(A^ / 2 - 1), Y.{N I 2)f (38) 

Y_,. = [Y.rnas,,(\\Y^a,,,{2),...Jn„a,,.{NI2-\)f (39) 

5 Comparing equations (36) and (12), in equation (36) we have to perform 

2iV' real-valued multiplications, and in equation (12), we need to perform IN x (N/I) 
complex-valued multiplications. Since one complex-valued multiplication requires four 
real-valued multiplications, the total number of real-valued multiplications is reduced by a 
factor of two in equation (36) compared to equation (12). 

10 A third, non-limiting, example embodiment of an MBAEC is shown in 

Fig. 6. The Hi matrix of size {NxN) is replaced by a vector Vi 94 of size (N x 1). A 
matrix also includes single column or row vectors, although die term vector is used here to 
distinguish from the N xJV matrix of the embodiments described above. Furthermore, 
the matrix Wi is replaced by another matrix Z,- 96 of size (JV x iV) . In this example 

15 embodiment, the previously transmitted frequency domain symbol Xm is multipHed m 

multipher 90 by a complex exponential term to compensate for the cyclic prefix, which has 
a length of L samples, i.e., e'^'^"' . The product output firom multiplier 90 is subtracted 
from Xi, and the difference is multiplied by matrix Zi 96. The vector and matrix products 
are summed in summer 98 to estimate the echo. 

20 The following equations describe the algorithm: 

Y,=V,.X,+Z,(x,-X,_,)t (40) 

where 

[x,-:l=[X,_,Le^- ,forA = 0,...,iV-l (41) 
X as the /:th transmitted frequency domain symbol, {N x 1) vector, defined above. 
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is the r.th estimated echo signal in frequency domain^ (N x 1) vector. 

Yz and is frequency domain vector and matrix of size (N x 1) and (N xN). 

Adjusting column vector and matrix with the LMS algorithm yields: 

V,,3=V,+//E,.(x:f (42) 
Z,,,=Z,+//E,(x:-X:J (43) 

where — Y? - is the /:th received frequency domain error, (N x 1) vector,// is the 
step size of the LMS algorithm, and • denotes element-by-element vector multiplication. 

Example procedures for the operation of the thitd example embodiment 
MBAEC shown in Fig. 6 are now described in conjunction with the routine 100 shown in 
flowchart form in Fig. 1. The signal to be transmitted is detected in the frequency 
domain (block 120). That signal is multiplied by the column vector Yi. The previously 
transmitted signal X^-/ is multiplied by a cyclic prefix compensating complex factor, which 
in this example embodiment is e^^"^^^^ , That product is subtracted from the current signal 
to be transmitted X^, and the difference is multiplied by the (N xN) matrix Zi (block 104). 
The products •X^ and Z^(x^ are summed to estimate the echo X O^lock 106), 

The estimated echo Y^ is subtracted from the actually received symbol Y^ to determine an 

error Ez (block 108). The coefficients of the column vector Yt and the matrix Z? are 
adjusted with the current error using the LMS algorithm (block 110). 

As a fourth example, non-limiting embodiment, the operations described in 
the second and third example embodiments may be combined. Specifically, the real and 
imaginary parts of the current and prior symbols are determined and used to simphfy the 
column vector and the matrix multiplication with matrhc in accordance with 

equations (23)-(39) set forth above. As a result, the computational complexity for 
performing the necessary data manipulations is reduced appropriately by a factor of four 
compared to the first example embodiment. 
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In a DMT application, the vectors X, , Y, , , and E, all have Hermitian 
symmetry. This property can be used to reduce the number of computations needed by 
only calculating the first N/2 + / coefficients of the vectors. Furthermore, the matrices 
and Z, need only be calculated and stored for the first N/2 + / rows. Moreover, the 
computational complexity may be further reduced by partitioning the vectors , X^_^ , , 
and into real and imaginary parts in accordance with equations (23)-(39). 

An MBAEC may be used m an asymmetric transceiver, such as an ADSL 
transceiver, with the modifications now described. When the modem's transmitter has a 
lower sampling rate than its receiver, the echo signal is mterpolated at the receiver (a fifth, 
non-limiting, example embodiment). When the modem's transmitter has a higher samplmg 
rate than the receiver, the echo signal will be decimated at the receiver (a sixth, non- 
limiting, example embodiment). 

For an interpolated MBAEC, assume the size of the IDFT in the transmitter 
is N and the si^e of the DFT in the receiver is ^N, where 8 m ADSL . Thus, the 
received symbol contains q times more samples than the transmitted symbol The received 
echo is therefore interpolated q times in the receiver. As one example interpolated 
asymmetric MBAEC, equation (40) is modified as follows: 

%^\.l^rep^^l.Xx^-%^,\ (44) 

where 

[X.. 1 - [X.„, ], - " . for ^ = 0. . . iV - 1 (45) 
Xi is the /:th transmitted frequency domam symbol, (N x 1) vector. 

lLrep^ IS equal to the vector Xi replicated q times, forming a block vector of 
contiguous vectors. For example, if ^ = ^ then is ^rept = [ X,^Xi^ X^^^. 

X is the fth estimated echo signal in frequenq^ domain, {qN x 1) vector. 

V, and Z, is frequency domain vector and matrix of size (^JV x 1) and 

{qN xN). 
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Adjusting and with the proposed LMS-algorithm yields: 

V,,,-V,+//E,.(xrei?:)' (46) 
Z„=Z,+^E.(x:-X:_J (47) 

where = Y. - Y, is the z:th received frequency domain error, ^ is the stepsize of the 
LMS algorithm, and • denotes element-by-element vector multiplication. 

Hence, the only differences from equation (40) are the replication of X, and 
the size of the matrices and the vectors. Some of the coefficients in the vectors and in the 
matrices will be small and can be neglected to reduce die computational complexity and 
the required memory. In addition, the computational complexity in the interpolated 
MBAEC may be further reduced by partitioning the vectors X, , X^,^^ , , 7, , and 
into real and imaginary parts in accordance with equations (23)-(39). 

For a decimated MBAEC, assume the size of the IDFT in the transmitter is 
qN and the size of the DFT in the receiver is N, where q- 8 'va. ADSL. Thus, the 
transmitted symbol contains q times more samples than the received symbol. This means 
that the received echo will be decimated q times in the receiver. Equation (40) is modified 
as foUows. 

Y, - Decimate^ ^ • X J, + Z, (x, - X,_, } (48) 



where 



2?r, 



[x._l-K,I.'^'\forA = O.....A^-l (49) 
X, is the /:th transmitted frequency domain symbol, {N x 1) vector. 

Decimate {(V,- • X ),q} means that the vector V, • X is decimated q times in 
the frequency domain. In other words the frequency domain vector V, • X. is blocked into 
q contiguous vectors of AT points each. These smaller vectors are then added together 
forming one vector of size (Nx1). 

is the 2:th estimated echo signal m frequency domain, (N x 1) vector. 
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Yi and 2i is frequency domain vector and matrix of size {^N x 1) and 

{NxqN). 

Adjusting Yi and with the proposed LMS-algorithm yields: 

V„,=X+//Ere/.,.(x:f (50) 
X,,,=Z,+juErep,{^:-XUj (51) 

where = - is the /:th received firequency domain error, ^repz is of size {qNx1) and 

is equal to the vector E^repHcated q times, forming a block vector of q contiguous vectors, 
ju is the stepsize of the LMS algorithm, and • denotes element-by-element vector 
multiplication . 

Again, the only differences to equation (40) are the decimation of • X? , 
the replication E^, and the size of the matrices and the vectors. Some of the coefficients in 
the vectors and in the matrices will be very small and can be neglected to reduce the 
computational complexity and the required memory. 

An asymmetric MBAEC routine (block 120) is now described in Fig. 8. A 
decision is made in block 122 whether the modem transmitter sampling rate is less than 
the sampling rate at the modem transceiver. If so, the symbol to be transmitted is 
interpolated at the modem receiver, e.g., at the analog-to-digital converter, and the 
estimated echo is calculated using the following equation: 

Y,=V,.Xr6p,+Z,(x,-X,,,) (52) 

(block 128). If the modem transmitter sampUng rate is greater than that for the modem 
receiver, the symbol to be transmitted is decimated at the receiver, e.g., at the analog-to- 
digital converter, and the estimated echo F is calculated in accordance with the following 
formula: 

Y, =Decimate{(V, .X,X^}+Z,(x, -X,_,) (53) 
For both the interpolation and decimation cases, the estimated echo Y^ is subtracted from 
the received signal to determine the error (block 132). The coefficients of the 
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column vector V, and the matrix are adjusted with the error using for example the 
LMS algorithm (block 134). 

A seventh, non-limiting, example embodiment of an echo canceller in 
accordance with the present ravention may be used in an asynchronous transceiver, such as 
an ADSL transceiver at the Central Office, with the modifications now described. In an 
asynchronous transceiver, the transmitted and received frames or symbols are not aligned 
in time. Hence, there is a time offset between the transmitter and receiver symbols. In 
this seventh embodiment, the echo canceller is de-coupled from the receiver. 
Nevertheless, asynchronous MBAEC may be applied to any one of the previously- 
described embodiments with a small modification. Instead of removing the echo in the 
frequency domain, the estimated echo signal is first transformed to the time domain by an 
IDFT operation, and then removed in the time domain sample-by-sample. The 
asynchronous MBAEC is illustrated in Fig 9. The estimated echo signal output from 
MBAEC is transformed mto the time domain by an TV - IDFT block 140 to produce an 
estimated echo in the time domain which is then removed in the time domain sample- 
by-sample firom the received echo . The time domain difference or error is 
transformed into the complex domain in the - DFT block which is located after the 
summation stage rather thaa before it as shown in Figs. 4 and 6. In addition, the 
computational complexity in the decimated MBAEC may be further reduced by 
partitionmg the vectors , X^_^ , 7^ , and into real and imaginary parts in accordance 
with equations (23)-(39). 

The proposed matrix-based adaptive echo canceller (MBAEC) has been 
implemented and simulated in different DMT-based transceivers such as a symmetric-rate 
DSL transceiver, an ADSL central office end transceiver, and an ADSL Subscriber end 
transceiver. Figs. 10 and 11 show the learnmg curve of the matrix-based adaptive echo 
canceller (MBAEC), implemented in a DMT-based Symmetric-rate DSL transceiver, using 
a simple lowpass model of the echo path channel The learning curve illustrates how the 
mean-square-error, MSE, is reduced during adaptation with the LMS-algorithm. The 
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simulation resiilts indicate that the present invention is reliable, accurate and performs 
according to theoretical expectations. 

The matrix-based adaptive echo canceller is weU-snited to follow slow 
channel fluctuations during transmission because the transmitted signal is used to adapt 
matrix coefficients. An echo canceller in accordance with the present invention is 
particularly efficient in a DMT system employing FDM (Frequency Division Multiplexing) 
because echo cancellation is only needed for a few carriers which leads to less 
computational complexity. Moreover, because the needed data are already available in the 
frequency domain, so there is no need to do any extra DFT or IDFT operations to 
implement the IVCBAEC, except for the asynchronous MBAEC, which requires one extra 
IDFT. Furthermore, some of the coefficients in the matrices and in the vectors will be 
zero or quite small, and therefore, can be ignored, which means that the computational 
complexity of and the required memory for the echo canceller will be reduced. 

While the present invention has been described with respect to a particular 
embodiment, those skilled in the art will recognize that the present invention is not limited 
to the specific embodiments described and illustrated herein. Different formats, 
embodiments, and adaptations besides those showQ and described as well as many 
modifications, variations, and equivalent arrangements may also be used to implement the 
invention. Therefore, while the present invention has been described in relation to its 
preferred embodiments, it is to be understood that this disclosure is only illustrative and 
exemplary of the present invention and is merely for the purposes of providing a full and 
enabling disclosure of the invention. Accordingly, it is intended that the invention be 
limited only by the scope of the claims appended hereto. 
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WHAT IS CLAIMED IS : 

1. An echo canceller for use in a transceiver configured to cancel a received 
echo signal completely in the frequency domain. 

2. The echo canceller in claim 1, comprising: 
a first matrix of coefficients^ and 

5 a second matrix of coefficients, 

wherein a combiaation of a product of the first matrix and a transmitted symbol 
and a product of the second matrix and a previously transmitted symbol is used to 
estimate an echo signal. 

3. The echo canceller m claim 2, wherein transmitted signals corresponding to 
10 the transmitted symbol and the previously-transmitted symbol are real-valued, and whereia 

the transmitted symbol and the previously-transmitted symbol are divided into real and 
itnagtaary parts before being combined with the respective matrix to reduce computational 
complexity. 

4. The echo canceller m claim 2, wherein the coefficients of the first matrix 
15 represent how an echo from a currently transmitted frequency domain signal affects a 

received signal. 

5. The echo canceller in claim 4, wherein the coefficients of the second matrix 
represent how an echo from a previously transmitted frequency domain signal affects the 
received signal. 

20 6. The echo canceller in claim 2, wherein the coefficients of the first matrix 

and the second matrix are adapted using a difference between the received signal and the 
estimated echo signal. 

7. The echo canceller in claim 6, wherein the coefficients are adapted using a 
least mean squares algorithm. 
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8. The echo canceller in claim 2, wherein the estimated echo signal is removed 
from the received signal. 

9. The echo canceller in claim 1, wherein the transceiver is a discrete multitone 
(DMT) transceiver. 

10. The echo canceller in claim 2, wherein the first and second matrices are N x 
N matrices, where N is a number of symbol samples. 

11. The echo canceller in claim 1, wherein a vector corresponding to a 
transmitted frequency domain symbol, a vector corresponding to a received frequency 
domain signal, and a vector corresponding to an estimate of the echo symbol are all 
Hermitian symmetric. 

12. The echo canceller in claim 1, comprising: 
a first vector of coefficients, and 

a second matrix of coefficients, 

wherein a combmation of a product of the first vector and a transmitted symbol 
and a product of the second matrix and a compensated, previously-transmitted symbol is 
used to estimate an echo signal. 

13. The echo canceller in claim 12, wherein the transmitted symbol and the 
previously-transmitted symbol are divided into real and imaginary parts before bemg 
combined respectively with the first vector and the first matrix to reduce computational 
complexity. 

14. The echo canceller in claim 12, wherein a compensation factor used to 
compensate the previously-transmitted signal is a complex exponential term. 

15. The echo canceller m claim 14, wherem the transceiver is a discrete 
multitone (DMT) transceiver and the compensation factor compensates for a cyclic prefix 
associated with the previously-transmitted signal. 
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16. The echo canceller in claim 1, wherein when a transmitter of the transceiver 
has a lower sampling rate than a receiver of the transceiver, the echo signal is interpolated 
at the receiver. 

17. The echo canceller in claim 1, wherein when a transmitter of the transceiver 
has a higher samplmg rate than a receiver of the transceiver, the echo signal is decimated at 
the receiver. 

18. An echo canceller for use in an asynchronous transceiver configured to 
cancel an echo signal, comprising: 

a first matrix of coefficients; 
a second matrix of coefficients; 

wherein a combination of a product of the first matrix and a currendy-transmitted 
symbol and a product of the second matrix and a previously-transmitted symbol is used to 

estimate an echo signal, and 

wherein the estimate of the echo signal is transformed to the time domain and then 

removed from a received signal in the time domaia. 

19. An echo canceller for use in an asynchronous transceiver configured to 
cancel an echo signal, comprising: 

a first vector of coefficients; 
a second matrix of coefficients, 

wherein a combmation of a product of the first vector and a currendy-transmitted 
symbol and a product of the second matrix and a compensated, previously-transmitted 
symbol is used to estimate an echo signal; and 

wherein the estimated echo signal is transformed to the time domain and dien 
removed from a received signal in the time domain. 

20. An echo canceller for use in a transceiver canceling an echo from a received 
signal in the frequency domain including circuitry configured to determme an estimate of 
the echo in the received signal using a frequency domain model of an echo path channel 
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that includes effects of interference and to subtract the echo estimate from the received 
signal. 

21. The echo canceller in claim 20, wherein the echo canceller is used in a 
discrete multitone (DMT) type transceiver and the interference includes inter symbol 
interference and inter-carrier interference. 

22. The echo canceller in claim 20, wherein the frequency domain model 
includes a first set of values that models how an echo from a currendy transmitted 
frequency domain symbol distorts the received signal and a second set of values that 
models how an echo from a previously transmitted frequency domain symbol distorts the 
received signal. 

23. The echo canceller in claim 22, wherein the first set of values is a first 
complex matrix and the second set of values is a second complex matrix, 

24. The echo canceller in claim 22, wherein the first set of values is a colimm 
vector and the second set of values is a matrix. 

25. The echo canceller in claim 24, wheretti the matrix is combined with a 
difference between the currently transmitted symbol and a product of the previously 
transmitted symbol and a compensating factor. 

26. The echo canceller in clakn 22, wherein the transmitted symbol and the 
previously transmitted symbol are divided into real and imaginary parts before being 
combined with the first and second sets of values, respectively. 

27. The echo canceller in claim 20, wherein when a transmitter of the 
transceiver has a lower sampling rate than a receiver of the transceiver, the echo signal is 
interpolated at the receiver. 

28. The echo canceller in claim 20, wherein when a transmitter of the 
transceiver has a higher sampling rate than a receiver of the transceiver, the echo signal is 
decimated at the receiver. 
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29. The echo canceller in claim 20^ wherein the transceiver is an as}Tachronous 
transceiver where transmitted and received frames or symbols are not aligned in time, and 

wherein the echo estimate is transformed into the time domain and then removed 
in the time domain from the received signal. 

30. A frequency domain echo canceller for use in a discrete multitone (DMT) 
type transceiver canceling an echo from a received signal in the frequency domain 
including circuitry configured to determine an estimate of the echo in the received signal 
using a frequency domain model of an echo path channel that includes effects of 
intersymbol interference and inter-carrier interference and to subtract the echo estimate 
from the received signal to provide a difference. 

31. The echo canceller in claim 30, wherein the frequency domain model 
includes a first set of values that models completely in the frequency domain how an echo 
from a currently transmitted frequency domain symbol distorts the received signal and a 
second set of values that models completely in the frequency domain how an echo from a 
previously transmitted frequency domain symbol distorts the received signal. 

32. The echo canceller in claim 31^ whereia transmitted signals corresponding to 
the currently and previously transmitted frequency domain symbols are real- valued. 

33. The echo canceller iti claim 32, wherein the currently transmitted symbol, 
the previously transmitted symbol, the received signal, and the difference are vectors 
having Hermitian symmetry. 

34. The echo canceller in claim 31, wherein the difference is used to adjust the 
first and second set of values. 

35. A method for reducitig an echo at a transceiver comprising: 

(a) combining ia the frequency domaia a currently transmitted signal with a first 
set of values resulting in a first combination; 

(b) combining in the frequency domain a previously transmitted symbol with a 
second set of values resulting in a second combination; 
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(c) combining the first and second combinations in the frequency domain to 
estimate the echo; and 

(d) using the estimated echo to reduce the echo in a signal received at the 
transceiver. 

36. The method in claim. 35, further comprising: 

determining a difference between the received signal and the estimated echo, and 
adjusting the first and second set of values using die difference. 

37. The method in claim 35, wherein the first set of values corresponds to a first 
matrix of coefficients and the second set of values corresponds to a second matrix of 
coefficients. 

38. The method in claim 35, wherein the first set of values corresponds to a 
column vector of coefficients and the second set of values corresponds to a matrix of 
coefficients. 

39. The method in claim 38, wherein the combining step (b) includes: 
multiplying the previously transmitted symbol by a compensation factor to produce 

a product; 

subtracting the product from the currendy transmitted symbol; and 
combining a resiilt of the subtractiag with the matrix. 

40. The method in claim 35, wherein when a transmitter of the transceiver has a 
lower sampling rate than a receiver of the transceiver, the method further comprising: 

interpolating the echo signal. 

41. The method in claim 35, wherein when a transmitter of the transceiver has a 
higher sampling rate than a receiver of the transceiver, the method further comprising: 

decimating the echo signal. 

42. The method in claim 35, wherein a vector corresponding to a transmitted 
frequency domain symbol, a vector corresponding to a received frequency domain signal, a 
vector corresponding to an estimate of the echo symbol are Hermitian symmetric. 
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43. The method in claim 35, further comprising: 

dividing the currendy transmitted symbol and the previously-transmitted symbol 
into real and knaginar}^ parts before the combining steps (a) and (b), respectively, to reduce 
computational complexity. 

44, The method according to claim 35, wherein the transceiver is an 
asynchronous transceiver, the method further comprising: 

transforming the estimated echo into the time domain, and 
removing the time domain, estimated echo signal from the received signal on a 
sample-by-sample basis. 
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ABSTRACT OF THE DISCLOSURE 

An echo canceller reduces an echo signal produced when the transmitted 
signal leaks back into the receiver via a hybrid. The echo canceller estimates the echo 
signal from the transmitted signal, and then subtracting the estimated echo signal from the 
received signal- In practice, the echo path channel m a DMT-modem is much longer than 
the cyclic prefix, and therefore, the received echo signal will be subjected to both ISI 
(inter- symbol-interference) and I CI (inter-carrier-interference). A traditional echo 
canceller, designed for a xDSL-modem, uses either a time domain adaptive FIR-filter or a 
combined echo canceller implemented in both time and frequency domain. A matrix- 
based adaptive echo canceller is implemented in the frequency domain. Various example 
embodiments are disclosed. 



Transmitter 



Transmitted signal 



r 



14 



Receiver 




Twisted-pair subscriber line 



Fi0. 1 



Copy the last L samples to beginning 




DI\/1T symbol (A/ samples) 



Cyclic prefix 
(L samples) 



CM 



CD 
-i— • 
- i — ' 

E 

CO 
05 



CO 
CO 



CO 
CO 







nsm 




CO 


[Z 







I 



CO 

VJ 



CM 
CO 

VJ 



CO 



D/A 




ii 


Add 


CP 


i 




P/S 



oo 

CM 



CD 
"O 
O 
O 

cz 
LU 




Encoder 



60 
VJ 



I 



N-IDFT 



Add 
CP 



Delay 



D/A 






1 r 



Matrix 


62 


Matrix 


Hi 




Wi 



Decoder 




64 



W-1 



Echo 
Path 

Channel 



N-DFT 




Remove 
CP 




A/D 



Pis. 4 



Fi£f. 5 ( MBAEC y -yo 



72 



Detect in the frequency domain a symbol 
to be transmitted (X/) 



I 



74 



Multiply X/ by matrix Hj and X/.^ by matrix Wj 



I 



/I 



76 



Add the products HjXj + WjXj.i to estimate the echo Yj 



1 



r 



78 



Subtract the estimated echo V^- from the actually received 
symbol to determine an error £/ 



I 



78 



Adjust the matrices W/^^ and VZ/^^ using the error £/, 
e.g., using the LI\4S algorithm 



Encoder 



60 
VJ 



I 



N-IDFT 




Add 


— ► 


D/A 




CP 





Delay 




(Nx1) vector 


94 




(NxN) Matrix 









96 




Z/[X/-X/.i(k)e'2^'-*/'^] 



98 



Echo 
Path 

Channel i 



Decoder 


- ^ "(5)- 5 


N-DFT 




Remove 




A/D 






CP 





Fi0. 6 



MBAEC ^ — ^100 



Detect in the frequency domain a symbol to be transmitted (X/) 



I 



•104 



Multiply Xf by Vector V/. Multiply X/., by a complex factor (e j^^^'^^'^), 
subtract the product from X;, and multiply the difference by matrix Z/ 



I 



^106 



Add the products X,V/ + Z\ [Xi - Xi.i(k) e'^^^^^^] to estimate the echo K/ 

i 



108 



Subtract the estimated echo V/ from the actually received 
symbol Y/ to determine an error £/ 



I 



110 



Adjust Vf and Z/ using the error £/, e.g., using the LMS algorithm 
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